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The operation of an interferometer for gravitational waves detection requires so- 
phisticated feedback controls in many parts of the apparatus. The aim of this 
lecture is to introduce the types of problems to be faced in this line of research. 
The attention is focused on the inertial damping of the test mass suspension of 
the VIRGO interferometer (the supcrattcnuator) : it is a multidimensional local 
control aimed to reduce the residual motion of the suspended mirror associated 
to the normal modes of the suspension. Its performance is very important for the 
locking of the interferometer. 

1 Introduction 

Operating an interferometer for gravitational waves detection requires the im- 
plementation of many active controls on the different parts of the apparatus. 
For instance, feedback controls are needed to reduce the laser frequency and 
power fluctuations, to keep the optical cavities in resonance and to maintain 
the interferometer output on a dark fringe. The basic idea is that the appa- 
ratus works at its best strictly around a well defined working point □ (laser 
fluctuations and shot noise at a minimum, optical cavities in resonance and 
interferometer output on a dark fringe). Internal and external disturbances 
overload the dynamic range of the interferometer. Therefore, the system has 
to be forced to remain in the correct working position. This is achieved via 
feedback controls. 

In this lecture we examine in detail one of the controls needed to reach 
the final goal of operating the Jnterferometer: the inertial damping of the 
VIRGO superattenuator (SA) a. The SA can be described as a chain of 
mechanical "filters", each one acting as a "spring" in 6 degrees of freedom. 
The normal mode frequencies of the SA range between 0.04 and about 2 Hz. 
At frequencies / >> 2 Hz the SA acts as a steep filter of the seismic vibrations 
of the ground (an attenuation factor of lO^^QlO Hz is expected). Therefore, 



♦LECTURE GIVEN AT THE INTERNATIONAL SUMMER SCHOOL ON EXPERI- 
MENTAL PHYSICS OF GRAVITATIONAL WAVES- URBINO (ITALY), SEPTEMBER 
6-18, 1999 



paper: submitted to World Scientific on January, 7, 2000 1 



in the interferometer detection band (10 Hz-few kHz), the suspended mirror 
is "disconnected" from the ground. Beside being a low pass filter for ground 
vibrations the SA is a tool for actively controlling the mirror position. Forces 
that move and steer the mirror can be exerted in 3 points of the S A chain: the 
inverted pendulum (IP) a, the marionetta (a special mechanical tool designed 
to steer the mirror) and the mirror itself (from a suspended reference mass). 

Keeping the interferometer in the operating position requires the mirror 
to have a maximum RMS relative motion of less the 10~^^ m. In the frequency 
range where the SA is fully effective (/ >> 4 Hz) the residual motion of the 
mirror is negligible. On the other hand, the mirror free motion in the region 
of the SA normal modes is ~ 100 fim . Feedback forces acting on the SA must 
reduce the mirror motion from 100 /xm to 10~^^ m. The dynamic range of a 
feedback system able to perform this control Jaas to be huge: the control is 
performed in three steps (hierarchical control □) . The first step is a damping 
of the SA normal modes in order to reduce the mirror residual motion to 
less than 10 /xm . This is necessary in order to control the mirror position 
acting on the lower stages without reinjecting noise into the detection band. 
We describe in the following an implementation of a high gain and wideband 
damping of the SA resonances. 

2 Experimental setup 

The setup (fig. ||) of the experiment is composed of a full scale superattenu- 
ator, provided with 3 accelerometers (placed on the top of the IP), 3 LVDT 
position sensors (measuring the relative motion of the IP with respect to an ex- 
ternal frame), 3 coil-magnet actuators. The accelerometers work in the range 
DC-400 Hz_and have acceleration spectral sensitivity ~ 10~^ ms~^Hz~^/^ 
below 3 Hz The sensors and actuators are all placed in pin-wheel configura- 
tion. The sensor and actuator signals are computer controlled by a ADC (16 
bit)-DSP-DAC (20 bit) system. The DSP handles the signals of all the sensors 
and actuators. It can combine them by means of matrices, create complex 
feedback filters (like the one of fig. ||) with high precision poles/zeroes place- 
ment and perform all the calculations at a high sampling rate (10 kHz). The 
suspended mirror is also provided with LVDT position sensors to measure its 
displacement with respect to ground. 

3 The approach to the problem of control 

The IP, where all the sensors are placed and on which the forces are exerted, 
has 3 main resonant modes: two translation {X, Y) and the rotation around 
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Figure 1. LEFT: the superattenuator; RIGHT TOP: logical scheme of the setup for the 
local active control; RIGHT BOTTOM: simplified view of the IP top table, provided with 
the 3 accelerometers. One LVDT position sensors and one coil-magnet actuator are also 
shown. 



the vertical axis (0). Each sensor is sensitive to all the modes and each 
actuator can excite all the modes. In control theory language, such a system 
is defined MIMO (multiple in- multiple out). Controlling a MIMO system can 
be very difficult. Our approach has been different: the signals of the 3 sensors 
are digitally mixed (using proper transformation matrices) to build up virtual 
sensors, sensitive to one mode only and "blind" to the others. At the same 
time, we build up virtual actuators able to excite each mode separately. The 
system is thus uncoupled into 3 SISO (single in-single out) subsystems. In 
terms of analytical mechanics this means the system is described in the normal 
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Figure 2. The logic of the diagonalisation: the output ui of the sensors are linearly combined 
by a matrix S in order to produce 3 virtual sensors outputs (xi), sensitive to pure modes. 
Three independent feedback filters Zi{s) are designed for the pure modes and 3 generalized 
forces qi are produced. The qi are turned into real currents {vi) to feed the actuators via 
the matrix D. 



modes basis. The equations of motions take the form: 

Xk +i^lxk ^Qk , 1...3 (1) 

where Xk is the fc— th normal mode, a;fc/27r is the corresponding resonant 
frequency and qk the generalized force on that mode. Let u — (-ui, ...,^3) the 
vector made by the outputs of the 3 sensors and let x — [xi, ■■■,xz) the vector 
made by the 3 virtual sensors. Analogously, let v the vector made by the 3 
currents driving the actuators and q the vector of the 3 generalized forces. 
The transformation from the system of the physical sensors/actuators to that 
of the virtual ones is operated by two matrices, such that: 

V = Dq (2) 
X = Su (3) 

The sensing (S) and driving (D) matrices are experimentally measured. The 
measurement procedure is described in details in ref. 9. The result of the 
measurement is determined by the mechanical characteristics of the system 
(resonant frequencies, quality factors), the geometry of the sensors and ac- 
tuators and their calibration, but it is not needed to know them in order to 
measure the matrices. The logic of diagonalisation is explained by fig. ^ 

After the diagonalisation the system can be considered as composed of 
3 uncoupled oscillators, and the control strategy for each of them has to be 
defined independently (see fig. ||). 
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Figure 3. Effect of the digital diagonalisation. LEFT: the output of the 3 virtual accelerom- 
eters when is excited. RIGHT: the output of the virtual accelerometers X and are 
compared. Different feedback strategies are required in the two cases, because X senses all 
the translational modes of the SA chain. 



4 Inertial damping: principle 

The control we describe here is called inertial damping because it is performed 
by using (mostly) inertial sensors (accelerometers). In the following , with 
the help of a simple model, we explain why this is the best choice to achieve 
a high performance damping. 

Let us consider a simple pendulum of mass m and length I. Let x be the 
abscissa of the suspended mass, 2:0 that of the suspension point. Let Ftb the 
external force on the pendulum (i.e. the feedback force to control it). The 
equation of motion is then: 

Ffh = rnx + 7i + k{x — xq) (4) 

where 7 is the viscous dissipation factor and k — mg/l. The control loop of 
such a system is sketched in fig. ^, where H{s) is the mechanical transfer 
function, G(s) is the compensator and out is the output of the sensor used. 
The goal of the control is to damp the pendulum resonance. This can be done 
easily with a viscous (theoretical) feedback force: 

i^fb = -7'i (5) 

Our sensors do not measure x. Their output is: 

^ ( X — xq for displacement sensors ,g, 
1 X for accelerometers 

Therefore, the actual "viscous" force that can be built if position sensors are 



paper: submitted to World Scientific on January, 7, 2000 5 




X 

Figure 4. The control scheme for a simple pendulum 



used has the form: 

K = -l^^i^-^o) (7) 

It can be easily shown that with such a feedback force the closed loop equation 
of motion (in Laplace space) reduces to: 

^(s) ^ ^ 9 7 77^ „ - ■ xo(s) (8) 

where Go = 7'/m is a gain parameter]^ measuring the intensity of the viscous 
feedback force, and Q is the open loop quality factor. When the loop is closed 
a damping of the resonance is achieved: 

Q' ""^^ 9 (9) 

Nevertheless, as the gain is increased, a larger amount of noise is reinjected 
off-resonance. This is associated to the term "Gos" in the numerator of 
and depends on the fact that the sensor used to build up the feedback force 
measures the position of the pendulum with respect to ground. Therefore, 
an infinitely efficient feedback would "freeze" the pendulum to ground (which 
is seismic noisy), reducing its motion at the resonance, with the drawback of 
bypassing its attenuation properties above resonance. 

The situation is fairly different when an inertial sensor is used. In this 
case the viscous feedback force is obtained by integrating the accelerometer 
output, and the output does not depend on xq: 

F,l = -7' / Mt (10) 



"^In a real feedback system a frequency dependent gain function G{s) rather than a gain 
parameter has to be considered. 
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Figure 5. Damping of a simple pendulum: the closed loop transfer function x(s)/xo{s) 
(magnitude) when a position sensor is used (LEFT) and when an accelerometer is used 
(RIGHT). 



The closed loop equation of motion is then: 

+ujq + (ljo/Q + Gojs 

A damping of the resonance is obtained (exactly as in the previous case) but 
without reinjection of off-resonance noise. In fig. ^ a simulation of the closed 
loop transfer function x{s) /xq (s) is shown in the two cases. Up to now we have 
considered a simple viscous damping. It is possible to increase the bandwidth 
of the control if the feedback force contains a term proportional to x (the 
double integral of the accelerometer signal). The result obtained in this case 
is shown in fig. || 



5 Control strategy 

In this section we extend the principles of the previous section and describe 
the strategy to control the SA. 

The basic idea of inertial damping is to use the accelerometer signal to 
build up the feedback force. As a matter of fact, an infinitely efhcient feedback 
using only the inertial sensor information, would null the acceleration of the 
pendulum, but it would not do anything if the pendulum moves at constant 
velocity. Such a control would be unstable with respect to drifts. Therefore, if 
the control band is to be extended down to DC, a position signal is necessary. 
Our solution was a merging of the two sensors: the virtual LVDT (position) 
and accelerometer signals are combined in such a way that the LVDT signal 
(l{s)) dominates below a chosen cross frequency /merge while the accelerometer 
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Figure 6. Inertial damping of a simple pen- 
dulum when a position feedback is imple- 
mented. 
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Figure 7. Merging of displacement and ac- 
celeration sensors (simulation for a simple 
pendulum). 



signal (a(s)) dominates above it (see fig. |^ and ref. 7). The feedback force 
has the form|^: 

Ff^^G{s)[a{s) + d{s)] (12) 

where G(s) is the digital filter transfer function (see fig. ^) and e is the 
parameter whose value determines /merge- We have chosen /merge ~ 10 mHz 
(corresponding to e ~ 5 • 10~'^). This approach stabilizes the system with 
respect to low frequency drifts at the cost of reinjecting a fraction e of the 
seismic noise via the feedback. 

We describe in the following the feedback design for the 3 d.o.f., starting 
from the the translational ones. The virtual X and Y sensors show many 
resonant peaks (the modes of a chain of pendulums) and this requires a more 
sophisticated feedback strategy. The digital filter used to control the transla- 
tion modes (G(s)) is shown in fig. || (LEFT). It shows three main features: 

• for 0.01 < / < 2 Hz the gain is proportional to /^^. This corresponds 
to the case of fig. the accelerometer signal is integrated twice and the 
feedback force is proportional to x; 

• for / > 2 Hz the gain is proportional to f~^. The accelerometer signal 
is integrated once: the feedback force is proportional to the velocity and 
a viscous damping is achieved; 

''Actually, the LVDT signal l{s) is properly filtered in order to preserve feedback stability at 
the crossover frequency and in order to reduce the amount of reinjected noise at / > /merge- 
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Figure 8. LEFT: Digital filter used for the inertial damping of a translation mode {X). The 
filter slope is in the range 10 mHz< / <3 Hz, /"^ for / >3 Hz. The unity gain is at 4 
Hz. The peaks in the digital filter are necessary to compensate the dips in the mechanical 
transfer function (see the transfer function of the X mode in fig. H). RIGHT: open loop 
gain function (measured). The phase margin at the unity gain frequency is about 25°. 



• the peaks visible in the filter are necessary to compensate the correspond- 
ing dips in the mechanical transfer function (H{s)) of fig. H, in order to 
make the feedback stable. 

Fig. H (RIGHT) shows the open loop gain transfer function G{s)H{s). 

The damping strategy for the Q mode is simpler: the virtual sensor (fig. 
^, RIGHT) shows one resonance peak only and no dips: no compensation is 
necessary. Apart from this, the feedback strategy is similar to the ones used 
for the translational modes. 



6 Inertial damping: experimental results 

The result of the inertial control (on 3 d.o.f.) is shown in figure |^. The 
measurement has been performed in air. The noise on the top of the IP is 
reduced over a wide band (10 mHz - 4 Hz). A gain > 1000 is obtained at 
the main SA resonance (0.3 Hz). The RMS translational motion of the IP 

(calculated as a;RMs(/) = Jp i^(i^)di^) in 10 sec. is reduced from more 

than 30 to 0.3 /im. The closed loop floor noise corresponds to the fraction of 
seismic noise reinjected by using the position sensors for the DC control and 
can, in principle, be reduced by a steeper low pass filtering of the LVDT signal 
at / > /merge and by lowering /merge: both this solution have drawbacks and 
need a careful implementation. 
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Figure 9. Performance of the inertial control {X, Y, © loops closed) of the superattuenuator, 
measured on the top of the IP: the left plot shows the acceleration spectral density as 
measured by the virtual accelerometer X (translation). The right plot shows the effect of 
the feedback on the RMS residual motion of the IP as a function of the frequency. 



Preliminary measurements of the displacement of the mirror with respect 
to ground have been performed in air, using an LVDT position sensor. The 
residual RMS mirror motion in 10 sec. is^: 

2^RMs(0.1 Hz) < 3 ^m. (13) 

When the damping is on such a measurement can provide only an upper 
bound because the LVDT output is dominated by the seismic motion of the 
ground. 

7 Conclusions 

In this lecture we have tried to outline how to face the problem of reducing the 
free motion of the suspended optical components of the VIRGO interferome- 
ter, associated to the resonances of the suspension. This is only the beginning: 
once the motion of the mirrors is reduced to a few microns, the lower control 
stages can operate to lock the interferometer in the correct operation state. 
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'^This number has been obtained with a simpler feedback design, less aggressive than the 
one of fig. ^ the gain raised as 1/f (pure viscous damping force), the crossover frequency 
was 30 mHz and no compensation of the dips was needed. 
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